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The energetic characteristics of two-shell carbon nanoparticles ("onions") with different shapes 
of second shell are calculated. The barriers of relative rotation of shells are found to be surprisingly 
small; therefore, free relative rotation of shells can take place at room temperature. The intershell 
orientational melting of the nanoparticle C6o@C'24o is studied by molecular dynamics. The param- 
eters of Arrhenius formula for jump rotational intershell diffusion are calculated. The definition of 
S . orientational melting temperature is proposed as the temperature when the transition probability 

over barrier between equivalent potential minima is equal to 1/2. The temperature of orientational 
melting of the nanoparticle Cqo®C24o is about 60 K. 
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I. INTRODUCTION 



C/3 . 

^ ■ The discovery of fuUerenes Q and the elaboration on method of their production in arc discharge |j give rise to 

(\ interest in another carbon nanostructures produced in arc discharge, in particular, nanoparticles with shell structure 

I ^^- A set of works is devoted to studying their structure and energetics [||- [|3[. Nevertheless, attention has not 

r^ yet been given to thermodynamical properties of carbon nanoparticles with shell structure. 

"*li , The melting of a single cluster can differ essentially from phase transitions in macroscopic systems p4|- [p2J. 

• ' Particularly, the melting of a mesoscopic cluster with shell structure can manifest itself as a hierarchy of rearrangements 

fx ] with breaking intershell orientational order and then breaking shell structure and order in particles positions inside 

•'-H shells. For example, in 2D mesoscopic clusters with Coulomb [Q- ||l^, screened Coulomb |2^, logarithmic |£l| and 

f;^.^^! dipole 1^,^ interaction between particles, the orientational melting (breaking the orientational order between the 

(— I 1 shells) precedes melting inside the shells. Namely, the reorientations of shells (jump rotational diffusion) or (with 

^ i' increasing temperature) free relative rotation of shells take place before shell structure breaking. This phenomenon 

'~~', is referred to as orientational melting. Moreover, the study of free relative rotation of shells may be interesting for 

_i ' nanomechanics [?]. 

^ ' The van der Waals interaction between atoms of neighbor shells in carbon nanoparticles is considerably weaker 

y—i than chemical bonds between atoms inside the shell. So it is natural that these nanoparticles are possible candidates 

T-H for orientational melting M . The possibility of orientational melting of long two-shell carbon nanotube was discussed 

^D pa . The orientational melting in carbon nanotube bundle was also theoretically studied pffl . 

O^ In the present paper the zero temperature energetic characteristics of two-shell carbon nanoparticle C6o@C'24o 
are calculated. The values obtained for barriers of relative rotations of shells are small enough for free rotation of 
shells to take place at room temperature. The orientational melting of this nanoparticle is studied here by molecular 
dynamics technique. The definition of orientational melting temperature is proposed. The corresponding temperature 



II. SIMULATION DETAILS 



Y^ for nanoparticle Ceo ©(^240 is calculated 

>>: 

The following reasons have determined our choice of nanoparticle shells. The TEM images show that the inner 

• '-H . shell of carbon nanoparticle can have a size that is close to that of fuUerene Cgo P^,p8| . The fuUerene Ceo with Ih 

^\ ' symmetry is the smallest fuUerene without adjacent pentagons in its structure. FuUerenes smaller than Ceo can not 

'^ , be directly extracted by the use of any solvent from soot, obtained in arc discharge (see, e. g., Refc. p^ , pO[ |). To 

■ - - ' explain this fact it was proposed that atoms of fuUerenes which belong to two adjacent pentagons can have chemical 

bonds with neighbor fuUerenes in soot [|32| . For example, chemical bonds between all neighbor fuUerenes are present 

in solid C36 [p4| . Therefore we consider Ceo as the smallest inner shell where the absence of chemical bonds between 

shells is very probable (it is a necessary condition for existence of relative rotation of shells). The single and double 

bonds lengths of Ceo used are 1.391 A and 1.455 A, respectively Q. We accept the fuUerene C240 with Ih symmetry 
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as outer shell of nanoparticle. This model gives the distance between shells in agreement with experiment |28 being 
close to the distance between graphite planes. Besides, fuUerene C240 with Ih symmetry have greater binding energy 
than fullerenes C240 with other structures . Several sets of geometric parameters corresponding to different shapes 
of fuUerene C240 obtained by ab initio calculations of minima of binding energy (5[|g,|9|] are used. Different shell 
shapes B, C, D and E were found by optimization of all independent geometric parameters of fuUerene C240 with I^ 
symmetry. The B and D shapes corresponding to global and local minima found by York et al Q that are close to 
sphere and truncated icosahedron, respectively. Shape E corresponds to the single minimum found by Osawa |q|. It 
is intermediate between shapes B and D. Shape C is rather close to shape E. It corresponds to the minimum found 
by Scuceria . The shape A is obtained by optimization of fewer of independent geometric parameters so that all 
atoms of this shape are arranged on the sphere |^ . 

We describe the interaction between atoms of neighbor shells by Lennard- Jones potential U = 4e((cr/r)^^ ~ ia/r)^) 
with parameters e = 28 K and a — 3.4 A. These parameters were used for the simulation of solid Cgo [£6[- "^^^ 
interaction between atoms inside shells are described by Born potential: 

where u^, u^ are displacements of atoms from equilibrium positions, ry are distances between atoms. We take 
a = 1.14 • 10"^ N/m and (3 = 1.24 ■ 10^ N/m. Born potential with these values of force constants gives an adequate 
internal vibrational spectrum of Cgo [0 ■ Born potential is correct only near the bottom of potential well. Nevertheless 
we believe that this potential is adequate for our simulation because we use it at temperatures that are one-two order 
of magnitude less than the temperature of fuUerene destruction. 

We studied the orientational melting of nanoparticle C6o@C24o with shape D of C240 by molecular dynamics 
technique. The simulations are performed in microcanonical ensemble. The equations of motion were integrated 
using the leap frog algorithm. We used the integration step r = 6.1 • 10~^^ s (about one hundred steps for period 
of atoms vibration inside shells). Initially the system has been brought to the equilibrium during 300-500 ps that 
is about 30-50 librations of shells. The average fluctuations of the total energy and temperature of the system fall 
and flatten out during this period. Then the system was studied during 100 ps. The average fluctuations of the 
total energy of the system were within 0.3 % and the average fluctuations of temperature were within 1.3 %. The 
angular velocities of shells change rather slowly to average the properties of system over the different directions of 
angular velocities during one computer experiment. Therefore all investigated quantities were averaged over 34-46 
different realizations of the systems at the same temperature but with different random angular velocities of shells 
corresponding to their distribution at temperature studied. 

III. RESULTS AND DISCUSSION 
A. Ground state energetics 

The global and local minima of total nanoparticle energy are found by optimization of three angles of their relative 
orientation. The total nanoparticle energy includes the energy of interaction between shells and the energy of shell 
deformation. We describe the relative orientations corresponding to minima of total energy in terms of three angles 
az^ cty and ax of subsequent rotations of first shell around axes OZ, OY and OX of coordinate system. The centers 
of both shells coincide with the center of coordinate system. The angles az, ay and ax were measured from the 
initial orientation shown on Fig. 1. Due to the high I^ symmetry of shells the number of any equivalent minima 
(global or local) is 60. Such equivalent minima correspond to different relative orientations of shells. The energies 
of interaction between shells and angles of one of the orientations corresponding to global and local minima of total 
energy of nanoparticle are listed in Table 1. 

The energies of interaction between shells calculated here are slightly less than 16.9 ||l^, 18.57 |0 and 20.3 |l^] 
meV/atom obtained using another representations of van der Waals interaction and are about three times less than 
estimation 65.3 meV/atom for graphite p^. Note, that the energy of total interaction between shells is not maximal 
for perfect sphere in comparison with other shapes of C240 contrary to the assumption of Lu and Yang [ |ll| . 

We observed that the angles of orientations corresponding to global and local minima are determined by the shape 
of second shell. For shapes C, D and E of C240 the initial relative orientation of shells (where symmetry axis of shells 
coincide) corresponds to global minima of total nanoparticle energy (note, that all these shapes of C240 are close to 
the truncated icosahedron). Several global minima for shape D are shown on Fig. 2a. One type of local minima is 
found for these shapes of C240 . For the shape B (which is close to sphere) orientations with coinciding symmetry axes 



correspond only to local minima (see Fig 2b). No minima correspond to such orientations for the "spherical" shape 
A. For the " spherical" shape of C240 two types of local minima are found. The differences AEioc in total nanoparticle 
energies between global and local minima are very small and also determined by the shape of second shell (see Table 
2). The differences AEioc decrease with decreasing the average deviation < ARi2 >=< \Ri2— < Ri2 > | > of second 
shell from perfect sphere, where Ri2 is the distance between an atom of second shell and the center of nanoparticle. 
Thus the differences AEioc are small for the "spherical" shape A and close to sphere shape B. The differences 
AEioc also decrease when the average distance between shells h —< Ri2 > — < Rn > approaches the distance rmin 
corresponding to the minimum in pair interatomic potential. This fact can be explained as follows: the smaller is the 
difference between h and rmin the lesser part of distances di2 between two atoms of neighbor shells corresponds to 
steeply rising interatomic potential well. Consequently, the change of distances di2 during relative rotation of shells 
causes the less change of interaction energy between shells. 

The calculated energies of shell deformation are presented in Table 2. The influence of shell deformation on the 
barriers of relative rotation of shells is studied as an example for barriers B^ of shell rotation around fivefold axes. 
(Barriers B5 were calculated for the relative orientation where symmetry axes of shells have the same directions). 
Comparison of barriers B5 calculated with and without shell deformation gives a difference less than 1 % for all five 
shapes of C240 investigated here. (Note that the barrier B5 calculated here for the shape E of C240 is 12 % less 
than that obtained by Osawa |^ who used the tandem of molecular orbital and molecular mechanics calculations). 
Therefore, the shell deformations are disregarded here in calculation of barriers of relative rotation of shells, i.e. 
lengths of bonds and angles between bonds inside shells are supposed to be fixed during intershell rotation. Note 



that an opposite situation take place, e.g. for clusters with logarithmic interaction between particles |21|. In this 
case, the interparticle interactions inside shell and between shells are the same and, therefore, the considering of shells 
deformation is necessary in calculation of barriers for rotation. The relative displacement of the centers of symmetry 
of shells causes an increase in intershell interaction energy. Therefore, the common center of symmetry of both shells 
is also supposed to be fixed during rotation. 

The barriers of relative rotation of shells in the nanoparticles under consideration are calculated for relative orien- 
tations corresponding to global minima of total nanoparticle energies. It is found that the values of barriers obtained 
for rotation are surprisingly small (see Table 2). Moreover, these barriers are only several times greater than barriers 
Ba in dependencies of interaction energy between only one atom of the second shell and the whole first shell vs. 
angle of rotation. For example, for the nanoparticle with shape D of C240 the barrier for rotation around fivefold 
axis is 158.8 K. Simultaneously, the maximal barrier among the barriers Ba for different atoms of the second shell 
is 21.6 K. Detailed analysis shows that maxima of barriers Ba for individual atoms in the same shell correspond to 
different angles of rotation and so the dependence of total energy on angle of rotation is essentially smoothed (see 
Fig. 3). Magnitudes of barriers of relative shell rotation are very sensitive to the shape of C240 and decrease when 
< Ai?i2 >— >■ and h — > rmin (analogously to the differences AEioc in interaction energies between global and local 
minima) . Note that the using of spherical shape of C240 leads to significant underestimation of barriers for rotation. 

The radii of shells of nanoparticle C6o@C24o are very close to radii of shells of (5,5)@(10,10) two-shell carbon 
nanotube. It is interest that barriers for relative rotation of shells per one atom calculated here for all considered 
nanoparticles are order of magnitude less than appropriate barrier in (5,5)@(10,10) two-shell carbon nanotube calcu- 
lated by Kwon and Tomanek P5[ . 

B. Molecular dynamics simulation 

We have investigated by molecular dynamics technique the angular velocity autocorrelation function of shells, the 
spectrum of shell librations, the frequency of shell reorientations, distributions of Eiler angles of relative orientations 
of shells, heat capacity of nanoparticle and barriers in intershell interaction energy corresponding to shell reorientation 
events. 

The dependence of total energy on temperature is used to calculate the heat capacity of nanoparticle. In the 
30 — 150 K temperature region investigated the heat capacity per one degree of freedom has no difference from the 
heat capacity of harmonic oscillator system within the accuracy of calculation that is less than 5 %. Only three 
degrees of freedom are accounted for relative orientation of shells. Therefore, as was to be expected, there is not any 
peculiarities in the dependency of heat capacity on temperature and the orientational melting of two-shell carbon 
nanoparticle has a crossover behavior. 

The dependence of shells reorientation frequency v vs. temperature T is shown on Fig. 4. The jump orientational 
intershell diffusion takes place where kT <C Bcf, B^f is an effective energy barrier of reorientation. We interpolate 
the reorientation frequency v for jump orientational intershell diffusion at temperatures 30 — 100 K by the Arrhenius 
formula (thick line on Fig. 4): 



^ = f^oexp(^-|^), (2) 

where Jig is a frequency multiplier. The fitting by least square technique gives B^f — f67 ± 22 K and 
fig = 540 ± f80 ns~^. Using a shorter temperature range T = 30 — 75 K for interpolation is found to have only 
a slight influence on calculated parameters B^f and Qq. Note, that effective barrier of reorientation B^f is in good 
agreement (within the accuracy of calculation) with the minimal Bmin and average Bav barriers for rotation of shells 
at zero temperature. Therefore it is possible to use barriers Bmin and Bav as effective barrier of reorientation to 
estimate the temperature of orientational melting of carbon nanotubes and nanoparticles with shell structure. 

The exponential increase of reorientation frequency v ends at temperatures 100 — 150 and this shows the beginning 
of free rotation of shells. It can be shown that the reorientation frequency v at temperature kT ^ B^f can be 
estimated by the expression 



where n is an average number of reorientations over the period of relative shell rotation (n « 5), /i and I2 are moments 
of inertia of 1-st and 2-nd shells, respectively. The dependence of reorientation frequency on temperature defined by 
Eq. (H) is shown on Fig. 4 by thin line. 

The prominent smooth distributions of Eiler angles of relative orientations of shells (Fig. 5), the disappearance 
of maxima in the angular velocity autocorrelation function of shells (Fig. 6) and in the spectrum of shell librations 
(Fig. 7) confirm that the free rotation of shells determines the thermodynamical behavior of the nanoparticle at 
temperatures greater than 140 . 

The temperature dependence of the " experimental" barriers Bre in intershell interaction energy corresponding to 
shell reorientations events is shown on Fig. 8. The barriers B^e are averaged over all observed reorientation events at 
corresponding temperature (30-70 reorientation events for each temperature investigated from range T = 40 — 55 K 
and 200-600 reorientation events for each temperature investigated from range T = 70 — 150 K). At temperatures 30- 
100 K, where the jump orientational intershell diffusion takes place, barrier Bre is in agreement (within the accuracy 
of calculation) with the minimal barrier Bmin for rotation of shells at zero temperature and with effective barrier of 
reorientation B^f. At temperatures 100 — 150 , where free rotation of shells begins, the magnitude of "experimental" 
barrier Bj-e grows and the increase of this barrier SB^e runs to 50 K at temperature 154 K. The increase dBj,e is 
greater than the dispersion ABav = 10 K of barriers for rotation of shells at zero temperature. Consequently, the 
increase (5 -Bre of the barrier can not be explained by climbing over the barrier with increasing temperature not only 
at their lowest point. Therefore, we believe that the increase AB^e of the barrier is the result of shell deformation. 
The increase of the dispersion ABre of "experimental" barrier with increasing temperature (see Fig. 9) also indicates 
the influence of shell deformation on this barrier. Note that energy of shell deformation is three order of magnitude 
greater than the increase 5Bre of "experimental" barrier in the result shell deformation. (The energies calculated of 
shell deformation Ed are in agreement with virial theorem Ed — kT{iN — 6)/2, where N is number of atoms in a 
shell.) 

Phenomena of order breaking in the systems with finite number of particles occur as a rule at some temperature 
range. This leads to problems in attempting to define the temperature of order breaking (see, e.g., Ref. ^^l ^'^d 
references herein). In the majority of cases a melting of clusters occurs with change of cluster structure. That is 
the system become spend time with increasing temperature not only in the ground state but also in states with the 
structure corresponding to other minima in potential energy of system. A mesoscopic system can fluctuate between 
different states separated by an energy barrier: the solid state corresponding to the ground state of the system and 
the liquid- like state corresponding to other minima in potential energy (see, e.g., Refc. [?]). To characterize such 
a system at its melting the quantity K{T) — 7; /7s was introduced pfl|, where 7; and 7^ are the probabilities that 
the system is at the temperature T in liquid-like state and solid state, respectively. In this case the temperature Tc 
corresponding to K(Tc) = 1 may be considered as melting temperature. 

The situation is different for the considered nanoparticle C6o@C'24o- The reorientations of shells are transitions 
between states corresponding to equivalent minima of potential energy. That is the structure of system does not 
change during orientational melting. To characterizes the melting of such systems, where diffusion occurs during the 
melting but the structure of the system does not change, we introduce the quantity K{T) = vt/^t, where vt is the 
frequancy of transitions between equivalent minima of potential energy and ujt is the frequancy of such oscillation 
wherein movement of particles directed along the path of transition. For considered systems we propose to define the 
temperature T^ corresponding to K{Tc) = 1 as the melting temperature. In this case a half of appropriate oscillations 
follows by transition to equivalent minima. Note, that proposed definition corresponds to a temperature of short 



order breaking and has no anology with phase transitions in macroscopic systems (contrary to above definition for 
the systems with transitions between inequivalent minima). 

The appropriate quantity characterizing the orientational melting of nanoparticle C6o@C'24o is Ko{T) — v/lo, where 
V is the shell reorientation frequency and w is the frequency of relative librations of shells. The equality Kq ~ 1 implies 
that a half of relative librations of shells begins in one minimum in dependence of potential energy of nanoparticle on 
angles of relative shell orientation and ends in neighbor equvalent minima. The estimation with the help of simulation 
performed gives the temperature Tc of orientational melting for nanoparticle C6o@C'24o with shape D of the second 
shell Tc « 60 K [Tc corresponds to Ko{Tc) = 1). Here the frequency of relative librations of shells is determined from 
the maximum in the spectrum of shell librations (see Fig. 7) . 

As we have shown, the barriers for rotation are very sensitive to the shape of shells. Therefore, the realization of 
possible orientational melting in many-shell nanoparticles is determined by their shape. The nanoparticles obtained in 
arc discharge are faceted in shape |g|,Q. However, their shape changes to almost spherical one when they are subjected 
to very strong electron irradiation in a high-resolution electron microscope p8| , p8||4l] . Accurate ah initio calculation 
of geometric parameters of large shells is necessary for performance of theoretical studies of possible orientational 
melting of many-shell nanoparticles. Nevertheless, the theory does not provide accurate coordinates. Some works 
predict that many-shell nanoparticles are faceted R,M , some that they are spherical |ll|,p^ , and some the transition 



from faceted to spherical shape for shells containing more than 3500 atoms J 42 L The calculations have also shown that 
the faceted nanoparticles transform to spherical under high temperature [nftl2J. Therefore, the barriers for rotation 
may decrease with increasing of temperature due to change of shell structure. Thus it is found that temperature of 
orientational melting Tc « 60 K for two-shell carbon nanoparticle is at least one order of magnitude less than the 
temperature of total melting. Analogously orientational melting can occur also in many-shell nanoparticles and short 
many shell nanotubes [p3. 

The carbon nanoparticles with shell structure are not the single example of different types of atom interaction 
inside shell and between shells. A two-shell spherical nanoparticle from M0S2 was produced |l^. We believe that 
orientational melting can also take place in nanoparticles from this and analogous sandwich materials {MX2^ M = 
Mo,W,X = S,Se). 

The orientational melting in a single nanoparticle may be revealed by IR or Raman study of the temperature 
dependence of width of spectral lines. The last must have Arrhenius-like contribution in reorientational phase (anal- 
ogously to the behavior in plastic crystals, see, e.g., Ref. Es] and references herein). Moreover, this study can give 
the estimation of reorientational barriers. Besides, NMR line narrowing can be observed in reorientational phase. 
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Table 1. 

The energies Eint of interaction between nanoparticle shells and one of the relative orientations of shells corre- 
sponding to the global and local minima of total nanoparticle energy; a^, ay and a^ are the angles of subsequent 
rotations of inner shell from initial orientation around axes OZ, OY and OX, respectively. 



Shape 


Eint 7 

(meV/atom) (i 


ctz cty a.x 
n radians) (in radians) (in radians) 


A 
A 
A 


15.034 
15.033 
15.032 


0.0819 
0.2495 
0.6283 


0.1452 
0.8128 
0.4634 


0.0540 

-0.0081 

0.0 


B 
B 


15.124 
15.101 


0.6283 
0.0 


0.4634 
0.0 


0.0 
0.0 


C 
C 


15.180 
15.098 


0.0 
0.6283 


0.0 
0.4634 


0.0 
0.0 


D 
D 


13.819 
13.777 


0.0 
0.6283 


0.0 
0.4634 


0.0 
0.0 


E 
E 


15.166 
15.061 


0.0 
0.6283 


0.0 
0.4634 


0.0 
0.0 



Table 2. 

The characteristics of second sheU shape: the average deviation of second shell from perfect sphere < ARi2 > 
and the difference between average intershell distance h and the distance rmin corresponding to the minimum in 
pair interparticle potential I = h — rmin', the differences A.Eioc in total energies of nanoparticle between global and 
local minima; the minimal and average barriers for rotation Bmin, Bav ± A_Ba^,, respectively, where the barrier Bav 
is averaged over all directions of rotation axis and ABav is its dispersion; the average energies of shell deformation 
Edi ± AEdi and Ed2 ± AEd2 for first and second shells, respectively, where the energies Edi and Ed2 are averaged 
over all relative orientations of shells and AEdi and AEd2 are their dispersions. 



Shape 


< ARa > 


I 


AEioc 
(°K) 


(°K) 


Bav ± ABav 

(°K) 


Edi ± AEdi 
(°K) 


Ed2 ± AEd2 
(°K) 


A 


-0.245 


0.0 


3.2; 5.5 


19.0 


20.5 ± 0.8 


2.09 ± 0.02 


34.56 ± 0.12 


B 


-0.258 


0.057 


76.7 


82.9 


122.1 ± 12.1 


1.62 ± 0.07 


29.98 ± 0.50 


C 


-0.289 


0.152 


287.4 


349.3 


363.1 ± 8.8 


2.17 ± 0.26 


18.19 ± 0.42 


D 


-0.119 


0.244 


144.4 


160.3 


177.3 ± 9.6 


3.75 ± 0.20 


34.40 ± 0.55 


E 


-0.299 


0.147 


368.3 


441.2 459.9 ± 12.9 


4.58 ± 0.44 


13.78 ± 0.38 



Captions for illustrations. 

Fig. 1. The fragments of two shells (shape D of second shell) at their initial orientations. OX, OY and OZ are 
axes of coordinate system. One fivefold axis of each shell is aligned with the axis OZ. One of the closest to axis OZ 
atoms of first and second shells (shown by black circles) lie in plane OXZ. This fixes the orientation of axes OX and 
OY. 

Fig. 2. The dependencies of binding energies for interaction between shells of nanoparticle on their relative 
orientation. Uz and ay are the angles of subsequent rotations of inner shell from initial orientation around axes Z and 
Y respectively. The angle of rotation around axis X is fixed equal to zero, a) shape D of second shell; b) shape B of 
second shell; 

Fig. 3. Interaction energies between first shell of nanoparticle and groups of atoms of second shell with shape D 
vs. angle a^ of rotation of inner shell from initial orientation around axis Z. An each group include all atoms with 
the same dependencies of interaction energy Ea between this atom and the first shell on angle of rotation. The curves 
corresponding to all 25 groups of atoms with different dependencies Ea for individual atom are shown by thin lines 
(23 groups from 10 atoms and 2 groups from 5 atoms). The dependence of total interaction energy between shells on 
angle a^ is shown by bold line. All energies are measured from their minima. 

Fig. 4 The dependence of shells reorientation frequency v on temperature T in Kelvin degrees. The interpolation 
by the Arrhenius formula at kT < Bre is shown by thick line. The estimation at kT > Bre is shown by thin line. 

Fig. 5 The distributions of Eiler angles 6, ip and (j> of relative orientations of shells at temperatures 21 K, 36 K 
and 140 K are shown by dotted lines, thin lines and thick lines respectively; a) the distribution of angle 0; b) the 
distribution of angle 9; c) the distribution of angle -0. 

Fig. 6 The angular velocity of autocorrelation function of the first shell at temperatures 21 K, 36 K and 140 K 
are shown by dotted lines, thin lines and thick lines respectively. 

Fig. 7 The spectrum of shell librations at temperatures 21 K, 36 K and 140 K are shown by dotted lines, thin 
lines and thick lines respectively. 

Fig. 8 The dependence of the " experimental" barriers Bre in intershell interaction energy on temperature T. 

Fig. 9 The dependence of the dispersion ABre of barriers in intershell interaction energy on temperature T. 
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